4. CONCLUSIONS 



Between the much-studied pure CDM and HDM models span a family 
of mixed dark matter models which can be parameterized by < fij, < 1, or 
equivalently, < mj, < 93/i^ eV. We are able to rule out all scale-invariant flat 
CDM+HDM models with n„ > 0.2 and Hq = 50 km s"^ Mpc"^ by comparing 
the abundance of collapsed objects at high redshifts in observations and in high- 
resolution TV-body simulations. This leaves only a narrow parameter space open 
since very small Q^, would recover the troubled CDM model. Is Q^, ^ 0.2 large 
enough to alleviate the problems encountered in the pure CDM model at low 
redshifts? Our simulations show that the Q^, = 0.2 model can better match the 
observed pairwise velocity dispersions and the two-point correlation function 
than the pure CDM model. We are analyzing other statistical quantities and 
completing the simulation to z = to further scrutinize this model. 
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3. LOW-REDSHIFT RESULTS 



While the high-redshift objects discussed above provide an upper bound 
on the neutrino fraction Q^, in CDM+HDM models, Q^, is bounded from below 
by various z ^ observations. After all, Q^, = leads us back to the troubled 
pure CDM model. (All low-redshift results discussed below are for the Q^, = 0.2 
model. We did not continue the simulation past z = 1.5 for the failed Q^, = 0.3 
model.) 

Figure 3 shows the two-point correlation function and the pairwise velocity 
dispersions at z = 0.25 (for Qrms-PS = 17.6 /xiiT). These are the two statistical 
quantities commonly used for comparison of models and observations. These 
should be compared with, e.g., Figures 8 and 12 for the pure CDM model in 
Gelb &: Bertschinger (CDMII 1994) before their halo break-up procedures were 
applied. In general, the particle velocities have been reduced from 1000-1200 
km s-i for the CDM model to 500-800 km s'^ for the n„ = 0.2 CDM+HDM, 
and the halo velocities have dropped from 500-600 km s~-^ to 200-400 km s~-^, in 
much better agreement with observations (Davis &: Peebles 1983; see also Zurek 
et al. 1994). The correlation function for the M > 2.5 X lO^^M© halos match the 
observed values fairly closely, although both the amplitude and the correlation 
length are too small for the M > 2.5 X 10^^ Mq halos. We are continuing the 
simulation to z = 0. It remains to be seen if there will be enough evolution 
between z = 0.25 and z = to bring ^(r) into agreement while keeping a 
reasonably low. 

In Figure 3, the smaller correlation amplitude and pairwise velocity for 
the hot particles compared to the cold ones (dashed curves) at r ^ 0.5h~^ 
Mpc reflect the free-streaming effect of the neutrinos. On larger scales, the two 
components track each other as expected. 



1000 



1000 






MM| 1 




" / 


\ 
\ 

^X 








--^ 


l/ 1 1 1 / 1 III 


mmI , 






0.1 



{h 



1 

Mpc) 



Fig. 3. The two-point correlation ^(r) and the pairwise velocity dis- 
persion a at z = 0.25 for the n„ = 0.2 CDM+HDM model. The 
dashed curves in both panels are for the CDM (top) and HDM (bot- 
tom) simulation particles; the solid curves are for dark halos with 
masses above 2.5 x lO^^M© (top) and 2.5 x lO^^M© (bottom). The 
dotted line indicates the observed ^(r) = (r/5/i~"^Mpc)~-^'® at z ^ 0. 
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Turner (1991) challenge hierarchical structure formation models, recent numer- 
ical studies (Katz et al. 1994) indicate that even the low-amplitude as = 0.5 
CDM model is compatible with the known distribution of high-redshift quasars. 

Since structure forms later in CDM+HDM models, a better constraint 
is provided by the amount of gas associated with the damped Lycc systems at 
z ~ 3 (Ma &: Bertschinger 1994 and references therein). Lanzetta (1993) and 
Lanzetta, Wolfe, &: Turnshek (1994) estimated that at 3 < z < 3.5, a fraction 
(2.9 ±0.6) X IQ-^h-^ of the critical mass density is in the form of dense neutral 
gas in these systems. Figure 2 shows their measured value and the predictions 
from our simulations for the two CDM+HDM models. The normalization of the 
power spectra corresponds to Qrms-PS = 17.6 /xiiT (Bennett et al. 1994). For 
comparison, we also plot the predicted Qg in the Q^, = 0.2 model for the higher 
normalization Qrms-PS = 21 /iK obtained by Gorski et al. (1994). The increase 
of a factor ~ 1.7 at the low mass end is too small to save the Q^, = 0.3 model. 

We conclude from Figure 2 that only the CDM+HDM models with n„ ^ 
0.2 can produce gravitational collapses early enough to account for the Qg in- 
ferred from damped Lycc systems at z ~ 3. 
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Fig. 2. The fraction of the critical density in baryons (hydrogen and 
hehum) associated with damped Lycc systems at z = 3.2. The thick 
curves assume the hydrogen column density threshold of Lanzetta 
et al. (1993), where all halos (with no overdensity constraint) were 
used. The thick solid curves are for f2j,=0.2 (top) and 0.3 (bottom) 
with Qrms-PS = 17.6 /xiiT (bottom), and the thick dashed curve is 
for f2j,=0.2 with Qrms-PS = 21 /xif for comparison. The top two 
dotted curves show Qg predicted by the f2j,=0.2 model, assuming 
all hydrogen in Denmax halos with Sp/p > 200 (bottom) and > 50 
(top) contributes to Lycc absorption; the bottom two dotted curves 
are for f2j,=0.3. (All dotted curves assume Qrms-PS = 17.6 /xiiT.) 
The shaded band shows the observed Qg from Lanzetta et al. 



of mass resolution on the mass functions at redshifts z ^ 1 is small. Since high- 
redshift objects are important for constraining our models, we carried out an 
additional test on the effect of resolution at high redshifts. Figure 1 shows the 
cumulative number densities of dark halos at z = 3 in two b = 2.5 (or as = 0.4) 
CDM simulations with different particle masses: 2.3 X lO^M© vs. 2.9 X lO^M©. 
Both simulations were performed in 100^ Mpc^ comoving boxes with Ho=50 km 
s"^ Mpc-^ The lower-resolution one (CDM16 in Gelb &: Bertschinger 1994) 
has 144^ particles and a Plummer force softening of 65 kpc (for h = 0.5); the 
higher-resolution one (Jain &: Bertschinger 1994) has a force softening of 20 kpc 
and 8 times more particles (and thus samples 8 times more wavenumbers ini- 
tially). This low-amplitude CDM model produces a comparable number of halos 
as the n„ = 0.3 CDM+HDM model at z ~ 3. Figure 1 shows only < 25% in- 
crease in the number densities when the mass resolution is increased by a factor 
of 8 and the force resolution increased by a factor ~ 3. Although we have not 
proven convergence, this test lends support to the validity of our mass resolution 
and high-redshift results discussed below considering that the uncertainties in 
cosmological observations {e.g. the normalization of the power spectrum) are 
generally larger than 25% and that we are ruling out models only based on 
order-of-magnitude discrepancies. 




Fig. 1. The effect of resolution on the halo number densities at 
high redshift. The Denmax algorithm was used to identify halos; 
no overdensity constraint was applied. The output at z = 3 from 
two b = 2.5 CDM simulations is shown. The thin curve is for 144^ 
particles with particle mass 2.3 X 1O"M0 and the thick curve is for 
288^ particles with mass 2.9 X 10^ Mq. An increase in mass resolution 
by a factor of 8 and force resolution by a factor of 3 makes only small 
difference in the number densities of halos. 

2. HIGH-REDSHIFT CONSTRAINTS 

Any hierarchical model of structure formation must account for the amount 
of collapsed objects observed at high redshifts. The quasar space density at red- 
shifts up to ~ 5 has been used to constrain the CDM model (Efstathiou &: Rees 
1988; Turner 1991; Katz et al. 1994). Although the analytic predictions by 
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ABSTRACT 

We report the most recent results from high-resolution numerical simu- 
lations of structure formation in two flat cold+hot dark matter models with 
neutrino mass densities Q^, = 0.2 and 0.3. We find that structure forms too late 
in all CDM+HDM models with Q^, > 0.2 to account for the amount of dense 
neutral gas in high-redshift damped Lyman-cc systems. The Q^, = 0.2 model at 
On z ^ provides a better match to observations than the pure CDM model. 

O 

^ 1. TV-BODY SIMULATIONS 

0\ Cosmological models with cold and hot dark matter (CDM+HDM) are pa- 

^ rameterized by the neutrino mass density Q^,. We chose to study structure for- 

mation in two f2 = 1 CDM+HDM models: (1) Q„ = 0.3, ficdm = 0.65, fibaryon = 
^ 0.05; (2) fij, = 0.2,f2cdm = 0.75, fibaryon = 0.05. The corresponding neutrino 
masses (assumed to be in one species for simplicity) are 7 and 4.7 eV, respec- 
tively. The Hubble constant is taken to be 50 km s~-^ Mpc"-^ (or h = 0.5). The 
^ non-linear growth of structure in the Q^, = 0.3 model has been shown (Davis 

Q.^ et al. 1992; Klypin et al. 1993; Jing et al. 1994) to provide excellent fits to 

observations at z = 0. However, the worry that galactic systems form too late 
in this model prompted us to investigate model (2), which assumes smaller fij,, 
I causing less free-streaming suppression and earlier formation of galaxies. 

2 We performed two high-resolution particle-particle particle-mesh (P^M) 

TV-body simulations to study the non-linear gravitational collapse in models (l) 
^ and (2). The same initial phases were used for the two simulations to give 

a controlled study of the effect of neutrino masses on galaxy formation. The 
primordial density perturbations were assumed Gaussian and scale-invariant. 
The simulation box was a cube with 100 Mpc comoving sides; the Plummer 
force softening was 50 kpc comoving (for Ho=50 km s~-^ Mpc"-^). A total of 23 
million particles were used: 2.1 million for the cold component, and 10 times 
more for the hot component to sample the neutrino momenta at a given position 
in phase space. The particle masses in the fij, = 0.2 model were 2.6 X 1O"^°M0 
for the cold and 6.6 X lO^M© for the hot component. The dynamic range and 
the number of particles make these the largest cosmological simulations ever 
published. 

Potential galactic sites were identified with the Denmax algorithm de- 
scribed in Gelb Sz Bertschinger (1994) and Ma Sz Bertschinger (1994). In ad- 
dition, an overdensity constraint of Sp/p > 200 was imposed to eliminate the 
diffuse, nonvirialized dark halos. 

The effect of force resolution in TV-body simulations has been studied ex- 
tensively by Gelb &: Bertschinger (1994), who also demonstrated that the effect 



